arXiv:1507.08845vl [astro-ph.SR] 31M2015 


Draft version August 3, 2015 

Preprint typeset using DT^X style emulateapj v. 5/2/11 


THE ABSOLUTE AGE OF THE GLOBULAR CLUSTER M15 USING NEAR-INFRARED ADAPTIVE OPTICS 

IMAGES FROM PISCES/LBT. 

M. Monelli 1,2 , V. Testa 3 , G. Bono 4,3 , I. Ferraro 3 , G. Iannicola 3 , G. Fiorentino 5 , C. Arcidiacono 5,6 , D. 
Massari 5,8 , K. Boutsia 3 , R. Briguglio 6 , L. Busoni 6 , R. Carini 3 , L. Close 9 , G. Cresci 5 , S. Esposito 6 , L. Fini 6 , M. 
Fumana 10 , J.C. Guerra 11 , J. Hill 11 , C. Kulesa 9 , F. Mannucci 6 , D. McCarthy 9 , E. Pinna 6 , A. Puglisi 6 , F. 
Quiros-Pacheco 6 , R. Ragazzoni 12 , A. Riccardi 6 , A. Skemer 9 , M. Xompero 6 , 

Draft version August 3, 2015 

ABSTRACT 

We present deep near-infrared (NIR) J, K s photometry of the old, metal-poor Galactic globular 
cluster M 15 obtained with images collected with the LUCI1 and PISCES cameras available at the 
Large Binocular Telescope (LBT). We show how the use of First Light Adaptive Optics system coupled 
with the (FLAO) PISCES camera allows us to improve the limiting magnitude by ^2 mag in K s . By 
analyzing archival HST data, we demonstrate that the quality of the LBT/PISCES color magnitude 
diagram is fully comparable with analogous space-based data. The smaller field of view is balanced 
by the shorter exposure time required to reach a similar photometric limit. 

We investigated the absolute age of M 15 by means of two methods: i) by determining the age from 
the position of the main sequence turn-off; and ii) by the magnitude difference between the MSTO 
and the well-defined knee detected along the faint portion of the MS. We derive consistent values of 
the absolute age of M 15, that is 12.9T2.6 Gyr and 13.3T1.1 Gyr, respectively. 

Subject headings: globular clusters: general globular clusters: individuakM 15 techniques: photometric 


1. INTRODUCTION 

Cosmological results based on recent Cosmic mi¬ 
crowave background (CMB) experiments (Boomerang, 
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are affected by an intrinsic degeneracy in the estimate of 
cosmological parameters, e.g. the Hubble constant Ho- 
To overcome this problem either specific priors_ or the 
resu lts of different experiments are used (Bennett et al. 


(BAO 

Eisenstein et al. 

|2005j), on supernovae observa- 

tions (J 

Riess et al. 1998, 

201l| and on gravitational lens- 

ing (Suyu et al. 2010, 2 

0T3) opened the path to the era 

of precision cosmology. However, the quoted experiments 


2014). 

Recent evaluations of the H 0 based on CMB provide 
values ranging from 70.0T2.2 km s -1 Mpc -1 (WMAP9 
Hinsh^we^ to 67 .8T0.9 km s -1 Mpc -1 (Planck 


Collaboration et al. 2015). Similar values have also been 


obtained by BAG plus supernovae using the so-called in¬ 
verse distance ladder sugge sting a value of 68.6T2.2 km 


s 1 Mpc 1 (Cuesta et al. 2015). On the other hand 


resolved objects (Cepheids plus supernovae) provide Ho 
values rangi ng from 73T2 (random) ±4 (systematic) km 
s -1 Mpc -1 ( Freedman fc Madore||2010| to 73.8T2.4 km 
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201 1|) . Slightly l arger values 
of the Hubble constant were obtained by |Suyu et al. 
(2013) using gravitational lens time delays (8O.CU4 7 km 


s -i Mpc -1 , uniform Ho in flat ACDM). 

The above estimates of the Hubble constant indicate 
that there is some tension between the results based on 
CMB and BAO and those based on primary and sec¬ 
ondary distance indicators. This critical issue has been 
addressed in several recent pa pers suggest ing a difference 
that range from almost 2 a (Efstathion 2014) to more 
than 2.5cr ( Riess et al.||2011 ). The quoted uncertainties 
on the Hubble constant open the path to new physics 
concerning the nu mber of relativistic species and/or the 


mass of neutrinos (Dvorkin et al. 2014 Wyman et al 


2014). Moreover and even more importantly, the above 


range in Ho implies an uncertainty on the age of the uni¬ 
verse -to- of the order of 2 Gyr. This uncertainty has 
a substantial impact not only on galaxy formation and 
evolution, but also on the age of the most ancient stellar 
systems, i.e. the globular clusters (GCs). 

The absolute age of GCs can be independently esti- 
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mated using stellar astrophysics and it is affected by the¬ 
oretical, empirical and intrinsic uncertainties. 

Theoretical- Stellar evolutionary models adopted to 
construct cluster isochrones are affected by uncertain¬ 
ties in the input physics. In particular, in the adopted 
micro (opacity, equation of state, astrophysical screening 
factors) and in macro-physics (mixing length, mass loss, 
atomic diffusion radiative levitation, color-temperature 
transformations). The impact that the quoted ingredi¬ 
ents have on cluster isochrones have been discussed in 


More recently it has been suggested to use as a refer¬ 
ence point a well defined knee along the low-mass regime 


detected in several old (lj Cen, 

Pulone et al. 1998| 

M4, 

Pulone et al. 1999; |Milone et al. 

2014; Braga etaT. 2015 

NGC3201, Bono etaT 2010; 47 

Tuc,|Lagioia et al. 2014 

NGC 2808, Milone et al. 2012; M71, |Di Cecco et al. 

2015) and intermediate-age (Sarajedini et al.||2009b| An 

et al. 2009b |a| 

I stellar systems and in the Galactic bulge 


detail in the literature ( 

VandenBerg et al. 

2013; Pietrin- 

ferni et al. 2004, 2009 Cassisi 

20l4|). T1 

he typical un- 


Off (MSTO)- is of the order of 10%. Thus suggesting 
that theoretical uncertainties does not appear to be the 
dominant source in the error budget of the absolute age 
of GCs. 

Empirical- The main source of uncertainty in the ab¬ 
solute age estimate of globular clusters are the individual 
distances (A/x 0 ~0.1 mag in the true distance modulus 
implies an uncertainty of 1 Gyr in the absolute age). The 
age estimate is even more affected when the uncertainties 
in the reddening cor rectio n and in the reddening law are 
taken into account (Stetso n et al.| [2014). 


Importantly, the massive use of multi-object fiber spec¬ 
trographs provided the opportunity to construct a firm 
metallicity scale inclu ding a significant fraction of GGCs 
( Carretta et al.|200 9), thus reducing the uncertainties in 


tne iron andln the ^-element abundances. 

Intrinsic- Dating back to more than forty years ago, 
spectroscopic investigations brought forward a signifi¬ 
cant s tar-to-star va riation in C and in N among cluster 
stars (|Osborn|1971|). This evidence was sou ndly comple¬ 
mented by variation in Na, Al 


and/or optical NIR CMDs. 

Nevertheless, one of the most difficult observational 
problem in measuring stellar magnitudes and colors in 
GCs is that they are intrinsically crowded stellar systems, 
and therefore the photometry of their stars is strongly 
limited by poor weather conditions. In particular, bad 
seeing (larger than ^ 1") has the effect of severely limiting 
the identification and measurement of faint stars. This 
means a systematic increase in the limiting magnitudes 
and in the photometric accuracy when moving from the 
outskirts to the innermost cluster center. 

Twenty five years ago the advent of the Hubble Space 
Telescope (HST) started a new era, and the high spatial 
resolution provided by space images collected in opti¬ 
cal bands allowed us to resolve the core of GCs. Re¬ 
cently, a similar resolution is becoming possible from 
ground based observations using NIR cameras available 
on 10m class telescopes assisted by adaptive optics (AO) 
systems. This technology allows ground-based observa¬ 
tions reach the diffraction limit over a modest field of 
view (^ lxl'). High-resolution NIR images of GCs can 
have a relevant impact on current astrophysical problems 
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Mg-Al (Suntzeff & Smith (T99lj Gratton et al. 2012). ([Ferraro et al. 2009; |B< 
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clusters (Bedin et al.|2004; Piotto et al 


2007) 


However, detailed investigations concerning the different (Neichel et al 
stellar populations indicate a difference in age that is^in 


GeMS/GSAOI available at the GEMINI-Sou th telescope 
~2014b|a~ ' 


Rigaut et al. 2014). This new 


canonical GCs, on average sh orter than 1 Gyr (Ventura 
et al.|200l Cassisi et al.|20 08 ). The intrinsic uncertainty 


does not seem to be the main source of the error budget 
of the GCs absolute age. 

To overcome or to alleviate the quoted uncertainties, 
have been suggested different approaches mainly based 
on relative age esti mates, th e so-calle d vert ical and hor- 
izontal methods flMarfn-Franch et al.|2009 Dotter et al. 
|2011| |VandenBerg et al. 20131) . In this context the rel- 
ative age is estimated as a difference between the clock 


system uses both natural guide stars for the tip tilt cor¬ 
rection and five artificial stars to close the loop and it 
has been able to delivery uniform NIR im ages approach¬ 


i ng h s diffraction limit in Galactic bulge (Saracino et al. 


2015) and halo (Turri et al. 2015) GCs. 


Tn this context we have collected AO images of the 
GC M15 (NGC 7078), using the Fir st Light Adaptive 


Optics (FLAO, Esposito et al. 2010) system mounted 
on the Large Binocular T elescope (LBT). This c luster is 


located at ^10 Kp c (10. 4, Durrell fc Harris|1993l 9.9 kpc, 
McNamara et al.||2004| 10.4 kpc, |van den Bosch et al.| 


2006) and is affected by moderate interstellar extinction 


horizontal branch or a specific point along the red gi- 

(E(B-V)=0.08 mag, 

Sandage et al.|1981 

0.10, 

Schlafly 

ant branch. The key advantage of these methods is that 

& Finkbeiner 

2011 0 

.12, Schlegel et al. 

1998 

). Most 


they are independent of uncertainties on cluster distance 
and reddening. However, they rely on the assumption 
that the reference points are independent of cluster age 
and introd uce new theoretical uncertainties (conductive 
opacities, Cassisi et al. 2007; e xtra-d eep mixing along 
the RGB, |Denissenkov & Weiss 2004). It goes without 
saying that the transformation of relative ages in abso¬ 
lute ages using a calibrating GC introduces the typical 
uncertainties already discussed. 


interestingly, it is a mong the most meta l-poor Galactic 
GCs ([Fe/H] « —2.4, Kraft fc Ivans||2003 ), and therefore 


it possibly traces the oldest component of our Galaxy 
(see Table 1 for the parameters assumed in our analysis). 
Notably, despite multi ple populations have been proved 
to ex ist in this cluster (Monelli et al.||2013 |Piotto et al.| 
2015), so far there is no evidence of multiple turn-off or 
sub-giant branches as for NGC 1851 (|Milone et al.||2008) 
that could affect the age estimate. 
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Fig. 1. — SDSS image of M 15 with superimposed the fields of the 
different data sets collected for this project. Detailed analisys of 
the central PISCES field will be presented in a forthcoming paper. 


2. OPTICAL AND NIR DATA SETS 

The present work uses four different data sets from 
different imagers. In the following we summarize the 
main properties of each of them. A summary is given in 
Table 2. 


2.1. LBT/PISCES data 

PISCES is an near-IR imager covering a wavelength 
range 1-2.5 /im with an Hawaii 1024 pxx 1024 px HgCdTe 
array, installed at the front bent-Gregorian focus (Mc¬ 
Carthy et al. 2001) of the LBT. It critically samples a 


diffraction-limited PSF with a plate scale of 0.0193 /px. 
Observations were carried using PISCES together with 
the FLAO system mounted on the DX (right) telescope 
of LBT. M 15 was observed on October 14-15 2011 during 
the Science Verification Time for the FLAO system. This 
is a twin of 672-actu ators voice-coil based, contactless 
adaptive mirror (Sali nari et al.|[l994 Davies et al. 2010) 


controlled by the means of a pyramid (Ragazzoni 1996) 
wavefront sensor. The FLAO uses solely Natural Guide 
Star as reference and it retrieves high Strehl-ratio over a 
broad wavelength range, reaching peak performance on 
bright reference in the NIR (80% in the H band). Once at 
regime phase the two FLAP systems will fe ed the LUCI1 
and LUCI2 cameras ( Le Fevre et al.||2003 ). 


TABLE 1 

Basic parameters of M 15. 


(m-M)o [mag] 

15.14 

Harris (1996 ); Durrell & Harris (1993) 

E(B-V) [mag] 

0.08 


Sandage et al. 

1981 


[Fe/H] 

-2.4 


Kraft & Ivans 

2003 



Two fields were acquired with the FLAO+PISCES 
setup, one centered on the cluster core, the other ap¬ 
proximately 3 arcmin South-West of the cluster center. 
In both cases the selection criterion for the field was the 
presence of a suitable star for wavefront analysis in the 
field-of-view, with magnitude R = 12.6 mag andR =12.9 
mag for the central and outer field, respectively. In the 
current investigation we will focus on the external field. 


A preliminary photometric analysis of the central field to¬ 
gether with a detailed investigation of the variation of the 
PSF across the field of v iew has already been discussed by 


Fiorentino et al. (2014). A comprehensive analysis will 
be addressed in a forthcoming investigation. During the 
observations, weather conditions were photometric with 


good natural seeing conditions (0.65 -0.9 , as recorded 
by the DIMM). The AO allowed to reach a mean FWHM 


// // 

of 0.05 and of 0.06 in the J- and in the LC s -band, re¬ 
spectively, as measured on the images. The Strehl ratio 
on the quoted images reached 28% (J) and 60% ( K s ) 
consistently with the expected scaling vs. wavelength. 


2.2. LBT/LUCI1 data 


One pointing with the spectro-imager LUCI1 at LBT 
was collected for calibration purposes. The set of obser¬ 
vations was secured on June, 21-22 2012 in the J and K s 


filters, under good seeing conditions (^0.7 ). The data 
set was taken with the center of the cluster in the NE 
quadrant of the image in order to include both the central 
and the outer field observed with PISCES. We adopted 
this observing strategy to constrain possible systemat- 
ics in the absolute calibrations due to positional effects. 
An off-source set of images was also taken to perform 
median -sky subtraction and superflat construction (see 
Section |2.4|). 


2.3. HST optical and NIR data 

Complementary data sets will be used in the analy¬ 
sis of the LBT images. In particular, a series of nine 
images has been retrieved from the HST archive. They 
were collected with the WFC3 in the F160W passband. 
Furthermore, we will make use of F606W and F814W 
photometry of M 15, retrieved from the “ACS Surve y of 
Galac tic Globular Clusters” database (Sarajedini et al. 
2007). 


Figure [l] shows an SDSS image of M 15 with super¬ 
imposed the footprint of the adopted cameras: LUCI1 
(black square), PISCES (blue), ACS (magenta) and 
WFC3 (red). The LUCI1 field (4'x4') was selected 
such to include both the center of the cluster and the 
PISCES field, which is located 2.7' from the center. Note 
the small field-of-view covered by the PISCES camera 
(21"x21 // ). The HST/ACS field is centered on M15, but 
due to the different position angle it is not overlapped 
with the PISCES pointing, which on the other hand 
falls into the HST/WFC3 set. Figure [2] presents a com¬ 
parison of the LBT/LUCI1 (left), HST/WFC3 (center), 
and LBT/PISCES (right) region corresponding to the 
full field-of-view of the PISCES camera. Other than the 
impressive improvement when comparing data from the 
same telescope but without (left panel) and with (right 
panel) assistance from the adaptive optics, it is clear that 
PISCES provide the best spatial resolution also when 
compared to the WFC3 (central panel). It is worth not¬ 
ing that the three images are stacked medians, thus the 
total exposure time is different in the three cases. More¬ 
over, while the left and right panels show images in the 
1C S , we only had available WFC3 images in the F160W 
filter, which is close to the H band. To highlight even 
more the PISCES performances, Figure [3] shows a zoom 
of the PISCES (left) and WFC3 (right) stacked images 
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already shown in Figure |2j The green squares mark a 
region of 5" x5". The image discloses at first glance that 
the number of sources is similar in both cases, suggesting 
a similar limiting magnitude (it is worth recalling here 
the shorter total exposure time in both PISCES J and 
IT S , see Table 2). On the other hand, the contrast in the 
PISCES images is by far better, and one can easily see 
that elongated sources in the WFC3 field are well sep¬ 
arated in the PISCES image, such as those next to the 
top left corner of the green square. 


TABLE 2 
Observations log. 


Telescope 

Sensor 

filter 

Exposures 

[s] 

Total Time 

M 

LBT 

LUCI1 

J 

13x20 13x40 

780 



K 

26x40 

1040 


PISCES 

J 

20x30 

600 



K 

42x15 

630 

HST 

WFC3 

F160W 

3x200 6x250 

2100 


2.4. Data reduction and photometry 

The acquisition and the basic reduction have been 
performed following a homogeneous approach both for 
PISCES and LUCI1 data: raw images have been secured 
by dithering the telescope within a 100 pixel random pat¬ 
tern to ensure good removal of bad pixels. Single images 
have been dark-subtracted, flat-fielded, resampled to re¬ 
move geometrical distortions, and registered. For the K s 
filter, a superflat obtained with the off-source sky images 
have been obtained and applied to the images to improve 
the low-frequency flat-field removal. 

The LUCI1 and WFC3 data have been ind epen ¬ 
dently reduc ed following the prescriptions of Monelli 
et al. (2010), and using a standard procedure based 
W~the DT^ UPHOTIV/ALLSTA R/ALLFRAME suite of 
programs (Stetson;|1987[ |1994|). Individual PSFs have 
been modeled for each image, using semi-automatic rou¬ 
tines. The input list of stars for ALLFRAME was gener¬ 
ated registering and matching the individual catalogues 
from single images. 

The case of the PISCES data deserved particular at¬ 
tention, because adaptive optics may provide PSFs char¬ 
acterized by spatial variations across the field due to 
anisoplanatism. This is especially true for the shorter 
wavelength J-band images, while the PSF in the K s is 
typically more stable even at the largest distances from 
the guide stars, as it scales with the isoplanatic angle 
and progressively with the wavelength. Therefore, to 
perform the photometry on these images we adopted the 


ROMAFOT suite of progra ms (Buonanno et al. 1983 r 
Bnonanno fc Iannicola|l989 ). The PSF photometry with 


ROMAFOTTs more lengthy when compared with sim¬ 
ilar packages available in the literature. However, it is 
has the key advantage of a graphical interface that al¬ 
lows the user to improve the local deconvolution of stel¬ 
lar profiles. The latest version of the code (Ferraro et 
al. 2015, in preparation) has been optimized to perform 
accurate photometry of crowded stellar fields on images 
collected with AO systems. In particular, it takes into 
account the spatial variation of the PSF across the field 


of view and the variation of the asymmetric, egg-like 
shape of the PSF. A preliminary discussion of the nu¬ 
merical algorithms and of the approach adopted to deal 
with the quoted i mage s has already been presented in 
Fiorentino et al. ( |2014 ). In passing we note that pho¬ 


tometry performed using asymmetric PSFs on ground- 
based (LUCI1@LBT) and space (WFC3@HST) images 
gives magnitudes that are, within the errors, identical to 
those measured using other photometric packages. 

The final adopted photometry was obtained in two 
steps. First, ROMAFOT was run on the mean J and K s - 
band images in order to create the master list of objects. 
The final photometry was obtained reducing the individ¬ 
ual images, and averaging the derived magnitudes. A 
Moffat analytic function was adopted to model the PSF, 
with /3= 2.0, 2.5 and <7=2.70, 2.05 for the J and K s im¬ 
ages, respectively. 

The final photometric catalogues were calibrated into 
the 2MASS photometric system in two steps. First, the 
LUCI1 photometric catalogue was calibrated using ^200 
stars in common with 2MASS and covering the entire 
field of view. Then three dozen local standards of the 
LUCI1 catalogue were used to calibrate the photometric 
catalogue based on PISCES images. In this context it is 
worth mentioning the crucial role that faint local stan¬ 
dards play in the accurate calibration of images collected 
with AO systems. The use of NGS for the tip tilt correc¬ 
tion and the modest field of view of current AO systems 
implies the selection of crowded stellar fields. In these 
cluster regions the photometric quality of the 2MASS lo¬ 
cal standards is quite poor, moreover, even the faintest 
2MASS stars are in these regions saturated. These are 
the reasons why the calibration of NIR images collected 
with AO systems does require a double step in the cal¬ 
ibration using 4-8m class telescopes to improve the lim¬ 
iting magnitud e^an d the_photometric accuracy of local 
standards (jBono et al. 2010). 

Keeping m mind the above caveats, Figure [4] shows 
the approach we adopted to perform the absolute cal¬ 
ibration. From left to right, the three rows show the 
residuals of the calibration as a function of magnitude 
(top: J; bottom: IT S ), the X and the Y coordinate. Small 
dots indicate the 2MASS stars in common with LUCI1, 
while the open squares show the ^35 LUCI1 stars used 
as local standard to calibrate the PISCES catalogue. No 
apparent trends are visible, and the residuals have null 
mean. 


3. RESULTS 

Figure [5] shows NIR CMDs of the central regions of 
M 15 based on four different imagers. From left to right, 
the (iC s , J — K s ) CMD based on: a) LBT/LUCI1 cover¬ 
ing a field of view of 16 arcmin squared (^17,000 stars); 
b) the intersection of space FT60IU-band (HST/WFC3) 
and ground-based iC s -band (LBT/PISCES) photometry 
covering a field-of-view (FoV) of 0.19 arcmin squared 
(^380 stars); c) LBT/PISCES covering a FoV of 0.19 
arcmin squared (^450 stars). The jump in limiting mag¬ 
nitude is compelling when moving from seeing-limited 
(LBT/LUCI1, b)) to adaptive optics (LBT/PISCES, c)) 
NIR images collected with the same telescope. Ind eed, 
the use of PISCES c amera together with t he FLAO |Es-| 
posito et al.| (|2012|); |Riccardi et al.| (|2010 ) allows us to 
move the limiting magnitude in the K-band down to 
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Fig. 2.— The sky region covered by the PISCES camera (right panel) as seen also by LUCI1 (left) and the WFC3 (center). Note the 
change in spatial resolution, from 0.118 /px (seeing limited), 0.13 /px (from space), 0.026 /px. 
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Fig. 3.— A zoom on the same PISCES (left) and WFC3 (right) images of Figure [ 2 ] The green square highlights the same sky region of 
5 // x5 // . The comparison clearly shows the gain in spatial resolution of the PISCES aata. 


^22.5-23.0 mag. To our knowledge the deepest iC s -band 
photometry ever performed in a GC. 

Panel d) shows a direct comparison between the CMD 
of the stars in common with PISCES (large black dots) 
and LUCI1 (green pluses). The two overlap over a mag¬ 
nitude range of ~ 3.5 magnitudes, from K s ?^16.6 mag to 
K s ^20.1 mag. The main sequence based on PISCES im¬ 
ages is narrower than the LUCI1 one, and no systematics 
appear in the comparison (see also Figure [3. Interest¬ 
ingly enough, the sequence of green symbols plotted in 
panel a) discloses a smaller photometric dispersion than 
the bulk of the main sequence stars in the LUCI1 field. 
The difference is mainly due to the low crowding level of 
these cluster regions. Indeed, the 450 stars measured in 
the PISCES FoV imply a density of 0.67 stars arcsec -2 
sq., that is one star every ^2,330 px 2 . This context is 
quite different when compared with the central pointing, 
since in this region the typical stellar density is a factor 
of ~50 higher. 

The CMD based on LUCI1 images allows the identi¬ 
fication of the typical evolutionary features of a GC. It 
covers more than 10 magnitudes in the K s band, and 
ranges from the tip of the RGB (K s ^9.5 mag) down 
to ^2 mag below the main sequence turn-off (K s ~18, 
J — K s ^0.25 mag). Moreover, M 15 also shows a well- 
populated horizontal branch ranging from K s ^16.5 to 


K s ^14 mag with the slope typical of NIR CMDs. Fi¬ 
nally, the RGB bump appears clearly in the luminosity 
function of the RGB at magnitude K s = 13.00 ± 0.05 
mag, indicated by the arrow. 

In spite of the good quality of LUCI1 photometry, it 
is thanks to PISCES and to the FLAO system that we 
have been able to identify, for the first time, the MSK 
(K s ^21.5 mag) in a very metal-poor GC. The quality 
of the CMD based on LBT/PISCES is further supported 
by the comparison with the CMD based on both space 
and ground-based NIR images. Data plotted in panel d) 
of Fig. [5] show that the limiting magnitude in F160W is 
similar to the limiting magnitude in the J-band. How¬ 
ever, the intrinsic error at fixed magnitude seems larger 
in F160W than in J-band CMD, and indeed the MSK 
cannot be easily identified in the above CMD. In pass¬ 
ing we also note that the exposure time in F160W is 3.5 
times larger than the exposure time in the J-band. 

We are dealing with photometric catalogues that have 
been collected using different telescopes equipped with 
different imagers and different sets of filters. Their pos¬ 
sible systematics might affect the absolute age estimates. 
To constrain this effect, Figure [6] and Figure [7] show the 
comparison of selected cluster isochrones with our data. 
A glance at the data plotted in this figure clearly shows 
the advantage of using CMDs based on optical and NIR 
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photometric data. Indeed, the optical-NIR CMDs do 
cover a range in color that is at least a factor of two larger 
compared with the optical ones. From left to right, Fig¬ 
ure [6] shows: a) the optical CMD based on ACS images 
(F6UBIF, F606W-F8UW); b) the optical-NIR CMD 
based on ACS and on WFC3 images (F606IF, F606W- 
FT60IU); c) the optical-NIR CMD based on ACS and 
on LUCI1 J-band images (F606VF, F606W-J ); d) the 
optical-NIR CMD based on ACS and on LUCI1 iGband 
images (F606VF, F606W-K S ). The cluster isochrone 
was com puted by adopting the ev olutionary tracks pro¬ 
vided byjVandenBerg et al. ([2014 ). The isochrones were 
transformed into the observat ional plane by adopting the 
color-temperature relation by jCasagrande fc VandenBerg| 
( 2014| ). We adopted an iron abundance of !FF/Hj=-2.4 


isE 

den et al.||1997[ 2t)00) and primordial helium content of 


(Kraft & Ivans 2003) an a-enhancement of a=+0.4 (Sne- 


Y =0.257 and a cluster age of 13 Gyr (red line). Cur 
rent theoretical framework is full y co nsiste nt wi th the 
set of isochrones adopted by Bono et alT] ([2010). We 
performed a series of test to constrain the optimal true 
distance modulus and reddening that provide a good si¬ 
multaneous agreement between theory and observations 
in the four CMDs plotted in Fig. 5. We found, using 
the Cardelli et al. (1989) reddening law, that a cluster 
reddening of E(B-V)=0.08 mag t ogether wit h a distance 
modulus of (ra —M) 0 =15.14 mag (|Harris 1996[ ,2010 edi¬ 
tion) do provide a good agreement in the quoted optical- 
NIR CMDs. Th e cluster reddening agrees with the 
value suggested by |Sandage et al.] ( |1981 ), but marginally 


smaller than the more recent estimates by |Schlafly fc 


Finkbeiner (2011) (E(B-V)=0.10). The true distance 


modulus agrees, within the erro rs, with diffe rent esti¬ 
mates available in the lite rature |Durrell & Harris 


1993 


van den Bosch et aL|2 006). Note that the assumption o 


a larger extinction would imply, at fixed age, a system¬ 
atic drift of the isochrones toward redder colors. The new 
discrepancy could be alleviated by a decrease in cluster 
age, but younger isochrones are characterized by a slope 
of the subgiant branch (SGB) that is too steep when 
compared with the data. This difference becomes more 
evident in the optical planes where the SGB is remark¬ 
ably thin. The above evidence indicates that adopted 
distance and reddening are mainly constrained by the 
morphology of SGB and RGB. The former one playing 
a crucial role, since it is les s affected by uncertaint ies in 
the adopted mixing length ( Salaris fc Cassisi||2Q15 ). 

The CMDs plotted in panels c) and d) also suggest a 
good agreement between theory and observations. In¬ 
deed, only a marginal shift in color is present at the 
base of the RGB in the F606VF,F606LF - J CMD. The 
anonymous referee also noted a similar shift in the RGB 
region between the HB and the base of the RGB. In the 
quoted cases, the isochrone are once again marginally 
redder than the observed RG stars. The above empirical 
evidence brings forward two relevant points. 

i)- Theory and observations disclose an overall very 
good agreement over more than ten magnitudes. The 
agreement becomes even more compelling if we take ac¬ 
count of the fact that we are dealing with optical and 
NIR data collected with space and ground-based facil¬ 
ities. This finding also supports the adopted chemical 
composition, the bolometric corrections and the color- 
temperature transformations together with the adopted 


true distance modulus and cluster reddening. Similar re¬ 
sults ha ve also been obtained in the literature by different 


groups (Dotter et al. 2008 

Sarajedini et al. 2009a[ Bres- 

san et al. 2012; DelkOmoc 

larme et al. 2012; VandenBerg 

et al. 2014) thus further supporting the current accuracy 


era of stellar astrophysics. 


ii)- The difference in color at the level of 0.03 mag 
between theory and observations has a marginal impact 
on the methods we are using to estimate the absolute 
cluster age. Indeed, they rely either on absolute (MSTO) 
or on relative magnitudes (MSK). 

Fig [7] shows a comparison with isochrones for different 
assumptions on the a enhancement and the metal con¬ 
tent. Left panel presents two isochrones of 12 (red line) 
and 13.5 Gyr (green). The effect of age appears clearly 
in the MSTO region, but as expected it does not affect 
neither the RGB nor the low main sequence. The right 
panel present two isochrones of 13 Gyr, with [Fe/H]=- 
2.2, -2.4 (red and green line, respectively) for different 
age and metallicity assumptions. 

4. THE ABSOLUTE AGE OF M15 

To estimate the absolute age of M 15 we devise here a 
double approach. The first is based on the MSTO posi¬ 
tion, and the second on the magnitude difference between 
the MSTO and the main MSK Table 3 summarizes the 
observables and the cluster ages based on the two quoted 
methods. From left to right the columns give the adopted 
CMD, the imager, the apparent magnitude, the color in¬ 
dex of the MSTO (m msto. Cl msto), and of the MSK 
(mMSK, CImsk ) together with the two cluster ages. 

4.1. The TO method 

The position of the MSTO is determined in the ob¬ 
servational planes the bluest MS point of the ridge line. 
Similarly, we estimate the magnitude and color index of 
the MSTO for a set of four isochrones of fixed metal¬ 
licity ([Fe/H] = —2.4, a = +04, Y = 0.25) and age 
between 12.0 and 13.5 Gyr (see £p3) once rescaled for the 
proper distance and reddening. This allows us to deter¬ 
mine that, at least in this age range, a linear relation 
exists between the MSTO magnitude and the age. The 
slopes of the above relations were obtained with a lin¬ 
ear Least Squares fit. The coefficients are listed in Table 
4, together with the predicted magnitudes of the MSTO 
and of the MSK as a function of cluster age. These slopes 
are a measure of the sensitivity of the MSTO as an age 
indicator in the different bands. For example, we derive 
that the F606IU band is ~1.8 times more accurate than 
the K s band. 

The age corresponding to M15 is derived by inter¬ 
polating the previous relations assuming the observed 
MSTO magnitude. The error budget has to take into 
account various sources, including error on the photo¬ 
metric calibration, the MSTO magnitude, the reddening 
and distanc^J In the case of the present data set, the 
photometric error varies depending on the filter used, 
from ^0.011 in the case of the F814W filter to ~0.04 

14 The error in the age is totally dominated by the propagated er¬ 
ror in the distance es timate. We note that the distance value from 
|van den Bosch et al. | (|2006|) provides error a factor of 2 smaller, 
wJnctL would imply a reduction of ~30% on the error on the abso¬ 
lute age determined with the MSTO. 
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Fig. 4. — Residuals of the calibrations as a function of magnitude (Left), X coordinate Center , and Y coordinate Right top and bottom 
rows refer to the J and K s bands, respectively. Large dots show the more than 200 stars in common between the LUCI catalogue and 
2MASS. Big open squares compare the ~35 stars in common between LUCI1 and PISCES, at significantly fainter magnitude. Note that 
the range of X and Y are different because of the size of the two cameras. 


mag for the LUCI1 K s . Nevertheless, while the differen¬ 
tial reddening is almost negligible along the line of sigh t 
of M 15 (E(B-V)=0.08=b0.01 mag, |Sandage et al.||1981|), 
the dominant source o f error is the distance (±0.15 mag, 


Durrell & Harris] 1993). We derive an absolute age rang¬ 
ing from 12.8T2.0 (E'©14 W band) to 14.0T3.0 Gyr (K s ), 
with weighted mean value of 12.9T2.6 Gyr. 


4.2. The MSK method 

The second approach is based on the magni¬ 
tude difference between the MSTO and the MSK, 
AM(MSTO — MSK). This method is based on the fact 
the magnitude and color of MSK in the low-mass regime 
of the MS are, at fixed chemical composition, essentially 


independent of cluster age {Bono et al. 2010). The key 
advantage of the MSK is that it is caused by collision- 
induced absorption (CIA) opacities of both H 2 -H 2 and 
H 2 ~ He i n the s urface of cool dwarfs ( Borysow et al.|200 r 


Borysow 2002). The MSK is independent of cluster age 


and anchored in a region of the MS that is marginally 
affected by uncertainties in the treatment of the con vec¬ 


tive r egime that is nearly adiabatic (Saumon & Marley 


2008). Recent empirical evidence indicates that the error 


m the absolute age of GCs based on this method are on 
average a factor of two smaller than th e canonical ones 
( Bono et aL||2010 Di Cecco et al.||20l5 ). 

'Therefore, this method otters a powerful observable to 
constrain the cluster age either as a color or as a magni- 
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Fig. 5.— (i<T s , J — K s ) CMDs of M 15. Panel a) - Data from LBT/LUCI1. The horizontal arrow marks the position of the RGB bump. 
The green pluses show the stars in common with the PISCES photometry. Panel b) - LBT/PISCES CMD of the outer M 15 field. Panel 
c) - Mixed CMD with HST/WFC3 [F160W] together with LBT/LUCI1 [K a \ data. Panel d) - Comparison between the CMD of the stars 
in common between LUCI1 (green pluses sa as in panel a)) and PISCES (large black dots). 


tude difference between the bend and the cluster MSTO, 
using both optical and infrared filter combinations. In 
the case of the present data set, we determined the po¬ 
sition of both points using the ridge line of the cluste r 
in the different CMDs (see e.g. Di Cecco et al. 2015). 
Note that lack of a sizable sample of MSTO stars m 
the PISCES photometry is a direct consequence of the 
modest field of view covered by the camera and by the 
radial distance of the pointing. To constrain on a more 
quantitative basis this relevant issue we selected the same 
cluster region covered by PISCES in the LUCI1 photom¬ 


etry and we found that it only includes 35 stars and the 
bulk of them are located at magnitudes fainter than the 
MSTO (see panel d) of Figure [BJ. 

To overcome this issue, we merged the LUCI1 and 
PISCES CMDs, and derived a unique ridge line. In 
this way, the TO region is sampled by the large num¬ 
ber of stars in the LUCI1 photometry, while the MSK 
is sampled by the PISCES deep catalogue. In particu¬ 
lar, following Bono et al. (2010), we define the color and 
magnitude of the MSK at the maximum curvature point 
in the low part of the MS. The outcome is shown in Fig- 

















The age of M 15 


9 



b) HSTV[ACS + WFC3] c) HSTVACS + LUCI J 


d) HST/ACS + L UCI K = 


a) HST/ACS 
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Fig. 6.— Optical and NIR CMDS based on HST data with superimposed the same isochrone for the labelled parameters. 


ure[8j where the combined LUCI1 and PISCES CMD are 
shown, together with the cluster ridge line, the MSTO 
and MSK. The approach adopted to compu te the c l uster 
ridge line has been discussed in detail by Di Cecco et al. 
(2015). The corresponding theoretical points were esti¬ 
mated using the same approach on the same isochrones 
adopted to estimate the cluster absolute age with the 
canonical MSTO method. 

Similarly to the age determination based on 
the MSTO, we derive linear relations between 
AM(MSTO — MSK) and the age (see Table 4). In 
this case, the method is slightly more sensitive than 
adopting the absolute value of the MSTO, by ^ 15%. 


The age is therefore derived interpolating the empirical 
value derived from the CMDs, and is reported in the 
last columns of Table 3. The estimated age ranges from 
12.9T0.9 (F606W band) to 13.7T1.0 Gyr ( K s ), with 
mean value 13.34=1.1 Gyr. 

4.3. Comparison with literature values 

The values shown in the last two columns Table 3 dis¬ 
close a general good agreement between the two derived 
age values. The estimate derived with the (K s , J — K s ) 
CMD is marginally higher than the other ones, but still 
well within the error bars. Also, the age derived with the 
MSTO are 0.3 to 0.9 Gyr smaller than the corresponding 
value derived with the MSK approach. 
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Fig. 7.— Comparison with isochrones of different ages (Zeft) and 
metallicities (right). 


The age of M15 has be en subject to a large num- 


& 


ber of i nvestigation s (e.g, _ 

Weiss 1998j 2002| |lVlcNamara et al.||2U(k| |T?> 
al.l I2005D . We stress that a straight compF 


et 


Salaris et al. 19971 Salaris 
Angeii 


‘ge 

n 


comparison of 


straigl 

literature estimates is complicated by the different theo¬ 
retical scenarios adopted. However, the age derived here 
is in good agreement with re cent estimates availa ble in 


the literature (12.8T 0.6 Gyr, Marm-Franch e t al. 2009, 
12.75T0.25T1.5 Gyr, VandenBerg et al.||2013[ ). Notably, 
these estimates are based on a difference approach. The 
former analysis use relative age of a population of clus¬ 
ters, anchored to an absolute scale using clusters with 
well-determined distance (NGC6752 via the subdwarf- 
based method). The latter uses the absolute magnitude 
of the MSTO. Note that we are using the same set of 
optical data as|VandenBerg et al.| (|2013|), and indeed we 
independently obtain a perfectly consistent value of the 
age based on the MSTO luminosity, close to 12.8 Gyr. 


5. DISCUSSION 

5.1. Comparison with NGC3201 and metallicity 
dependence 

In the present work we studied NIR photometry of 
the very metal-poor cluster M15. Nonetheless, it is 
interesting to extend the analysis to a more metal- 
rich regime in order to explore the dependency of 
the AM(MSTO — MSK) method on both the age a nd 
the metallicity of the target system. [Bono et al. 
(|2010|) presented a similar analysis for the cluster 
IN GC 3201 ([Fe/H] ~ —1.5). By adopting the same 
theoretical framework applied here we derive analo¬ 
gous linear relations correlating the magnitude difference 
AM(MSTO — MSK) with age, assuming [Fe/H] = —1.5. 
The obtained derivatives (e.g.: 0.13T0.01 mag Gyr -1 
and 0.08T0.01 mag Gyr -1 for the F606W and the K s 
band, respectively) are very similar to those derived for 


more metal-poor isochrones suitable for M 15. This sup¬ 
ports that the AM(MSTO — MSK) diagnostic can be 
fruitfully used over a wide range of metallicities. 

Moreover, by comparing the absolute position of the 
MSK in different isochrones, we find that: i) at fixed 
metallicity, the magnitude of the MSK changes by at 
most 0.02 mag, for ages larger than 10.5 Gyr; ii) simi¬ 
larly, at fixed age the MSK moves by ^0.02 mag, when 
moving from [Fe/H] = —2.4 to [Fe/H] = —1.5. Interest¬ 
ingly, these values seem independent of the wavelength, 
at least in the spectral range covered by the F606W to 
K s the passbands. Overall, the position of the MSK 
seems to be a reliable anchor, marginally dependent on 
the age and the metallicity, at least in the ranges in¬ 
vestigated so far. However, a more systematic theoret¬ 
ical investigation is needed, in order to to constrain in 
detail the sensitivity of the MSK, and in turn of the 
AM (MSTO — MSK) method, over the full range of ages 
and metallicities typical of Galactic GCs. 


5.2. Cosmological implications 

There is mounting evidence that we are in the era of 
precision cosmology. Recent estimates of the Hubble con¬ 
stant suggest a precision of the order of 3% (74.3T2.1 


-l 


-l 


Mpc U Freedman et al.||2012 73.8T2.4 km 
Riess et al.^2011 ). Plain physical arguments sug- 


km s 
Mpc -1 

gest, that in a flat Universe, the age of the Universe -to~ 
is connected with the Hubble constant -Hq-, the matter 
density parameter -f l m - and with t he d ark energy den¬ 
sity -tl\- by the following relation (De kel et al.| [l997): 


to = [1 - (On - 0.7 D a )/ 5.8]/(1.3 x h 0 ) (1) 

where ho = Ho/100 km s -1 Mpc -1 is the current espan- 
sion rate of the universe and to is the age of the Universe 
today in units of 10 10 Gyr. 

Using the recent estimates of the cosmological param¬ 
eters provided by WMAP and Planck (U m =0.315lo ois 

n A =0.685lg;flg,~ 


Hinshaw et al. 2013 


Planck Collabora- 


tion et al.||2014| ) we found t 0 = 10.70±0.82 Gyi£3 The 


above cosmological age is within lcr of current stellar 
ages. However, it is on the short limit if we also take 
account of the time for structure formation (z~8, t< 1 
Gyr). A similar approach to constrain the age of the 
Universe today is to use together with the above cosmo¬ 
logical parameters also the estimate of Ho provided by 
Planck+WMAP (Hq =67.3±1.2 km s -1 Mpc -1 ). It is 
worth mentioning that Ho is a prior in the CMB solu¬ 
tion ranging from 20 to 100 km s -1 Mpc -1 . Note that in 
this context the precision on H 0 is of the order of 1.8%. 
The age of the Universe we found is to=11.75±0.21 Gyr, 
while the detailed inversion of the CMB maps provides 
£o=13.817±0.048 Gyr. The above estimates appear in 
quite good agreement with the current absolute ages of 
oldest GCs, and we note that a smaller age of 10.71T0.50 
Gyr is derived assuming the H 0 value from Riess et al. 

I20TT1 

However, we are facing the evidence that the uncertain¬ 
ties on cosmological ages are systematically smaller than 


15 Note that we did not assume an Eisten-deSitter cosmological 
model (Q m =l, 12 a =0) because the age of the Universe we obtain 
is systematically younger (to=8.6±0.28 Gyr) than suggested by 
recent “stellar” and “cosmological” estimates. 
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TABLE 3 

Age determinations 


CMD 


Cameras 

m M STO 

CImsto 

rriMSK 

CImsk 

tMSTO 

tMSTO —MSK 

F814VE, F606VE 

- F8UW 

ACS 

18.801T0.011 

0.490T0.016 

21.642T0.036 

0.829T0.052 

12.8T2.0 

13.3T0.6 

F160VE, F814VE 

- F160W 

WFC3+ACS 

18.027i0.040 

0.652T0.041 

20.296T0.051 

1.045T0.059 

12.54=2.7 

13.44=1.3 

F160VE, F606W 

- F160W 

WFC3+ACS 

18.027T0.040 

1.147T0.041 

20.377T0.051 

1.838T0.061 

12.6=42.7 

12.9=41.3 

K S ,J-K S 


PISCES 

18.010T0.043 

0.238T0.044 

20.160T0.290 

0.581T0.299 

14.0=43.1 

13.74=1.4 


TABLE 4 

Magnitude v s ages derivatives. 


Mag 

dMag(TO)/dt 

dMag(TO-MSK)/dt 

F606VE 

0.09=40.01 

0.11=40.01 

F814VE 

0.08=40.01 

0.09=40.01 

J 

0.06=40.02 

0.07=40.02 

F160VE 

0.06=40.01 

0.07T0.01 

K s 

0.05=40.02 

0.06=40.02 

Mag 

dCI(TO)/dt 

dCI(TO-MSK)/dt 

F606W-F8UW 

0.01=40.01 

0.02=40.01 

F606W-F160W 

0.03=40.01 

0.02=40.01 

J-K s 

0.01=40.01 

0.01=40.01 


1 Gyr for estimates based on Hq and become smaller than 
a few hundreds of Myrs for CMB determinations. On the 
other hand, the uncertainties affecting the absolute age 
of GCs ranges from 2 (MSTO) to 1 (MSK) Gyr. The dif¬ 
ference in precision between cosmological and stellar ages 
is going to become even m ore prominent, since the next 


gene ration of experiments (Riess et al. 2011 Freedman 
et al.|2012 ) plans to improve by a factor of two the preci- 


sion on Hq. Taken at face values the “stellar” estimates 
do not allow us to validate the evaluations provided by 
cosmology. It is clear that accurate absolute ages for a 
sizable sample of extremely metal-poor GCs can shed 
new lights on this long-standing and intriguing problem. 


6. CONCLUSIONS AND FUTURE REMARKS 

We have presented new NIR data of the metal-poor 
globular cluster M15 obtained with the LUCI1 and 
PISCES cameras available at the LBT, and comple¬ 
mented with archival optical and NIR HST data. The 
analysis of the data raised important points: i) ground- 
based adaptive optics camera can compete with the HST 
in terms of photometric depth and resolving power in 
moderately crowded stellar fields, with smaller invest¬ 
ment of telescope time; ii) Tests performed with the 
ROMAFOT package suggest that the data reduction of 
images from adaptive optics instrument requires the de¬ 
velopment of novel techniques to model the complicated 
PSF of these imagers. 

The MSK is an important feature in the CMD of 
GCs that can help revising the age of these fundamen¬ 
tal stellar systems. The obvious advantage of using the 
AM(MSTO — MSK) approach, being a differential mea¬ 
surement, is that it is not affected by the errors either 
in the distance or in the reddening. This is reflected in 
the significantly smaller error bars in the age determina¬ 
tions. The analysis presented in this paper reveals that 
different photometric bands provide different sensitivity 
to the method. In this sense, the sensitivity decreases 
for increasing wavelength from the optical to the near 
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Fig. 8.— (K, J - K) CMDs from LBT data (black: LUCI1; red: 
PISCES) with superimposed the derived ridge line. The big dots 
mark the MSTO and MSK points. 


infrared. 

Our data analysis allowed us, using the PISCES data, 
to measure the MSK along the main sequence of M 15. 
We use two diagnostics to estimate the absolute age of 
this cluster: the magnitude of the MSTO and the magni¬ 
tude difference AM(MSTO — MSK). The two methods 
provide consistent results, and a mean absolute age of 
12.9T2.6 Gyr and 13.0T1.1 Gyr, respectively. 

A systematic theoretical analysis of the dependence of 
the MSK magnitude over a wide range of ages and metal- 
licities is mandatory to firmly establish the uncertainties 
affecting the AM(MSTO — MSK) method. Nonetheless, 
our results suggest that using high quality, optical data¬ 
bases such as those based on existing HST data can pro¬ 
vide a fundamental starting point to globally revise the 
age of the GC systems. The PISCES data presented 
here soundly demonstrate the potential of ground-based 
NIR data using adaptive optics technology to obtain deep 
photometry in crowded stellar fields. Moreover, current 
ongoing observing facilities at the 10m class telescopes 
using either SCAO (FLAO at LBT) or MCAO (GEMS 
at Gemini South), are providing excellent d ata reaching 
the MSK in many Galactic GCs (NGC 1815, Turri et al.| 


2014, Turri et al. 2015). The same outcome applies to 
the near future facilities like ERIS at ESO VLT. Nev¬ 
ertheless, a giant leap forward is foreseeable when the 
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next generation of extremely large telescope, equipped 
with NIR detectors and AO systems (e.g. Deep et al. 


2011 Greggio et al. 2012 Schreiber et al. 


2014), - wm 


the 


uro- 


be available: the GMT-Giant Magellan Telescop ^ 
TMT-Thirty Meter Telescop j^H, an d the E-ELT- 
pean Extremely Large Telescop^j). 

Realistic predictions according to up-to-date instru¬ 
mental specifications for the E-ELT+Maory+Micado 
configuration suggest that the predicted limiting mag¬ 
nitude (K ~27.2 assumin g a crowding leve l expected 
for the core of a GC, e.g. Deep et al. 2011|), is signif¬ 
icantly fainter than the expected MSJx magnitude for 
any stellar system in the Milky Way Halo (within 100 
Kpc, MSK^24.5). This is true also when restricting to 
a relatively short integration time of 600 s, comparable 
with that of the PISCES data presented in this paper. 
Preliminary analysis suggests the error bar on the age 


will be smaller than 1 Gyr for the entire sample of GCs. 
This means that we will be able to successfully apply 
this method, for the first time, also to the nearby Local 
Group galaxies and their GCs. 
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